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ABSTRACT 



The one thing among many that must be prepared for every flight is the making 
of a flight plan. Pilots must use charts or graphs from the flight manual to compute the 
fuel flow that is essential to a flight plan. Since this requires many steps of 
interpolation to compute the specific conditions that cannot be read directly from flight 
manual, it is time consuming and increases the probability of making a mistake. This 
problem obstructs the execution of various mission changes and continuous sorties. 

A computer program for personal computer or hand-held calculator is developed 
to compute the desired fuel flow by modifing the equations for an 'IDEAL' airplane. 
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I. INTRODUCTION 



A. BACKGROUND 

A pilot must be furnished with flight performance information for his aircraft 
that will permit him to operate the aircraft efficiently at different airspeeds, different 
altitudes, different values of gross weight and different external loadings. 

This information is collated from a combination of theoretical determinations 
and flight test data and presented to the pilot in a flight handbook in either a tabular 
format or in a series of graphs, or in combinations of these methods. 

B. FUEL FLOW RATE 

Probably the most important parameter in optimal performance of an aircraft is 
the rate at which fuel is consumed. It is obvious that when the fuel is exhausted the 
aircraft will no longer fly, so the pilot is very interested in the rate at which the fuel is 
being consumed. But also of importance is the fact that as the fuel is burned the weight 
of the aircraft is decreased and the performance of the aircraft is changed. It is 
therefore quite important that the pilot be able to forecast the rate of fuel consumption 
for each portion of each flight. 

Fuel flow rate for a turbo-jet aircraft is a direct function of the thrust required 
for the specified flight condition, and in equilibrium flight thrust required is equal to 
the total drag of the aircraft. 

C. TOTAL DRAG 

The total drag of an aircraft in equilibrium flight may be subdivided into induced 
drag (the drag associated with the generation of lift) and parasite drag (the drag 
associated with the resistance of the fuselage to motion through the air). The induced 
drag is a function of the altitude of the aircraft (air density) and the weight of the 
aircraft, and the parasite drag is also a function of the altitude as well as being a 
function of the profile drag of the aircraft. The aircraft drag, which varies with the 
configuration and external loading (weapons, fuel tanks, et cetera) is usually expressed 
in terms of a Drag Index (DI), with the zero Drag Index usually being a measure of 
the drag of the aircraft in the 'clean' configuration. 
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D. THRUST 



Although it is the thrust of the engines that provides the performance for a 
turbo-jet aircraft, as long as there is sufficient thrust available to provide the required 
performance, the pilot is more interested for the flight planning in the fuel flow rate. 

E. FLIGHT PLANNING 

To plan a flight, the pilot must determine what will be the loading of the aircraft 
(gross weight) and what effect this loading will have on the drag (Drag Index). The 
pilot then enters the flight handbook to determine the fuel flow rate for the airspeed 
and altitude at which he desires to fly. 

If the gross weight. Drag Index and airspeed are exactly the same as values listed 
in the handbook tables or exactly the same as the valued on the graphical lines, the 
determination of the fuel flow rate is simply a matter of reading (correctly) the values 
from the information present in the handbook. 

If, however, any of the input parameters are not identical to the handbook data 
points, the pilot must interpolate between the values presented. If weight. Drag Index, 
airspeed and altitude are all off-presentation points, the pilot may be forced to do 
several cross interpolations. This may even involve interpolations between values listed 
on different pages. Not only is this a time consuming and difficult task, but the 
possibility for error is expanded with each interpolation. 
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II. APPROACH TO THE PROBLEM 



A. STATEMENT OF THE PROBLEM 

A method was sought that would permit the pilot to rapidly obtain the fuel flow 
rate for various flight combinations of gross weight. Drag Index and airspeed. And, 
regardless of the ease of access to the information and completeness of the data 
obtained, this method must provide fuel flow rate at the same degree of accuracy as the 
flight handbook. 

In addition, it was planned that the method of data retrieval should be as simple 
as possible with an ultimate goal of using a hand-held computer for the solutions. 

B. BASE LINE APPROACH 

The standard method of using the flight handbook is for the pilot to search and 
follow several lines and tables of the handbook to determine the required fuel flow rate. 
This requires that altitude, gross weight and Drag Index must be considered, but this 
method is too long, takes too much effort and is open to errors. 

The procedure can be simplified by computerizing the flight manual. Most 
attempts at this method have involved the development of curve-fitting routines that 
will- supply numerical modeling of the handbook data. Although the end use of these 
equations is quite rapid, the development of these routines is a long and arduous task. 
Inasmuch as each equation usually fits only a small part of the data, the computer 
programs for such a method are usually quite lengthy. 

The approach undertaken in this project was the use of the equations for an 
'IDEAL' aircraft to develop simplified equations that can be used in a short computer 
program. 

Three steps were undertaken in the development of these equations: 

a) The fuel flow rate data in the handbook is converted to thrust required. 

b) The thrust required equations are adjusted so that they represent the actual 
aircraft. 

c) The corrected thrust required equations are then reconverted to flow rate as the 
program output. 
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C. DETAILED METHODS 
1. Concept of Kj and K 2 

Determinaiions are made for equilibrium, unaccelerated conditions where drag 
equals thrust. 

The basic drag equations for an ideal aircraft are used [Ref l:p. 160] 

Cj^_ = ^ AR 

D = (l/2)p V2 S Cd 

All of the principal items of flight performance involve steady state flight 
conditions and equilibrium of the airplane. For the airplane to remain in steady level 
flight, equilibrium must be obtained by a lift equal to the airplane weight. And the 
airplane drag equals the thrust required to maintain steady level flight. 

The total drag of the airplane is the sum of the parasite and induced drag. 
Parasite drag is the sum of pressure and friction drag which is due to the basic 
configuration and, as defined, is independent of lift. The equation of parasite drag can 
be expressed as a function of squared velocity. 

Dq = (1/2)Cd^ P S v2 

Let Kj = (l/2)Cj)^ p S 

Do - K, V2 

Induced drag is the undesirable but unavoidable consequence of the 
development of lift. 

In the airfoil, the local lift vector is aligned perpendicular the the local relative 
wind and hence is inclined behind the vertical by the angle a . Consequently, there is a 
component of the local lift vector in the direction of Vqq ; that is, there is a drag 
created by the presence of downwash. 

This drag is defined as induced drag, denoted by Dj [Ref. l:pp. 152-155] 

D- = (l/2)p S V2 Cq. 

Cj). = / (^ e AR) 

Di = (Cl) 2 /(27t e AR)p S 
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= [(2W/P S v2)2 / (271 e AR)] p S 
= (4w2 p S v 2) / (p2 n e AR) 

= [(2w2)/(pS7reAR)](l/v2) 

If the whole item of right side except the squared velocity is changed to some 
constant, this equation can be expressed function of squared velocity. 

Let K 2 = (2w2) / (p S Jt e AR) 

D- = K2(1/v2) 

While the parasite drag predominates at high speed, induced drag 
predominates at low speed. Figure 2.1 illustrates the variation with speed of the 
induced, parasite, and total drag for a specific airplane configuration in steady level 
flight. 

The curve of drag or thrust required versus velocity shows the variation of 
induced, parasite, and total drag. Induced drag predominates at low speeds when the 
airplane is operated at maximum lift-drag ratio, (L/D)j^^j^ the total drag is at a 
minimum and the induced and parasite drags are equal. 

The effect of a change in weight on the thrust required is illustrated by Figure 
2.2. The changes in the drag curves with variations of airplane weight, configuration, 
and altitude furnish insight for the variation of range, endurance, climb performance, 
etc., with these same items. 

First, the primary effect of a weight change is a change in the induced drag at 
any given speed. Thus, the greatest changes in the curves of thrust and power required 
will take place in the range of low speed flight where the induced effects predominate. 
The changes in thrust and power required in the range of high speed flight are 
relatively slight because parasite effects predominate at high speed. The induced effects 
at high speed are relatively small and changes in these items produce a small effect on 
the total thrust required. [Ref. 2;p. 99] 

Second, Figure 2.3 illustrates the effect on the curve of thrust required of a 
change in the equivalent parasite area f of the configuration. Since parasite drag 
predominates in the region of high flight speed, a change in f will produce the greatest 
change in thrust and power required at high speed. Since parasite drag is relatively 
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Figure 2. 1 Airplane Thrust and Power Required. 
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EFFECT OF WEIGHT CHANGE 



THRUST 

REQUIRED 

(LBS) 




0 100 200 300 4(X3 500 

VELOCITY-KNOTS 



POWER 

REQUIRED 




0 100 200 300 400 500 

VELOCITY-KNOTS 



Figure 2.2 Effect of Weight on Thrust and Power Required. 
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small in the region of low speed flight, a change in f will produce relatively small 
changes in thrust required at low speeds. The principal effect of a change in equivalent 
parasite area of the configuration is to change the parasite drag at any given airspeed. 
[Ref 2:p. 101] 

Third, a change in altitude can produce significant changes in the curves of 
thrust required. The effect of altitude on these curves provide a great part of the 
explanation of the effect of altitude on range and endurance. [Ref 2:p. 101] 

Figure 2.4 illustrates the effect of a change in altitude in the curves of thurst 
required for a specific airplane configuration and gross weight. As long as 
compressibility effects are negligible, the principal effect of increased altitude on the 
curve of thrust required is that specific aerodynamic conditions occur at higher true 
airspeeds. 

As being investigated above, the total drag, same as thrust required, is the 
function of squared velocity with different constant. 

Total drag equation is 

D = Kj v2 + K 2 / 

where Ki = Cr\ p S 
* ^0 

K 2 = (2 W2) / (p S 7t e AR) 

The first step of this project is to determine the Kj and K 2 values to obtain 
the thrust required. It can be seen that the Kj is the function of altitude and drag 
index, K 2 is the function of gross weight and altitude. So if Kj and K 2 values can be 
computed with any flight condition, the thrust required will be solved at any airspeed. 

2. Determining the Kj 

To determine the Kj, use the minimum drag concept. 

At the minimum drag point, the airplane can be flown maximum endurance 
and the induced drag and parasite drags are the same. 

The first step for determine the K| is to investigate the thrust equation. The 
thrust equation can be established with fuel flow for each airspeed and some factor C, 
THRUST = FUEL FLOW x C 

Since the fuel flow of the turbojet powered airplane is proportional to the 
airspeed, the thrust required is proportional to the fuel flow. 
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Figure 2.3 Effect of Equivalent Parasite Area,f, Thrust and Power Required. 
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EFFECT OF ALTITUDE CHANGE 
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Figure 2.4 Effect of Altitude on Thrust and Power Required. 
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The factor C is the value of the military power devided by fuel flow at military 

thrust. 

C = MILITARY THRUST / FUEL FLOW ( at military thrust ) 

The military thrust and fuel flow information for any condition can be 
obtained from NATOPS FLIGHT MANUAL. With this factor C, the thrust required 
for each speed can be computed by fuel flow, came from NATOPS FLIGHT 
MANUAL, multiplies factor C. Table 1 and Figures 2.5 and 2.6 show the relationship 
and the result between thrust and velocity. 

TABLE 1 

THRUST REQUIRED AND FUEL FLOW VERSUS VELOCITY 



VELOCITY(KTAS) 

360 

400 

440 

480 

520 

560 

600 

MIL 



FUEL FLOW(LBS/H) THRUST REQUIRED(LBS) 



7,725 
9,014 
10,623 
12,462 
14,652 
17,377 
21 393 
25,680 



6,558 

7,652 

9,018 

10,579 

12438 

14752 

18,160 

21,800 



Example) 



Thus 



military thrust ; 21800 lb 
fuel flow (at mil) ; 25680 Ibs/hr 
Factor C = 21800/25680 
= 0.84891 

Thrust required = 0.84891 x Fuel Flow 



Maximum endurance speed at specific gross weight, altitude, and drag index 
can be figured out with the graph of maximum endurance mach number in NATOPS 
MANUAL Figure 2.7. This true mach number needs to be transfered to KTAS to 
substitute into the drag equation. 

Next step is to compute the minimum drag. 

Minimum drag can be computed from the fuel flow at maximum endurance 
multiplied by the factor C. The maximum endurance fuel flow can be obtained from 
the NATOPS MANUAL Figure 2.8. 

And so the minimum drag is the twice of the parasite drag. 



D 



total 



= 2Dj = 



2Do = 



■mm 
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l igurc 2.C riirust Required for Velocity. 
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Figure 2.7 Maximum Endurance Mach Number. 
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Figure 2.8 Maximum Endurance Fuel Flow. 
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Thus the drag equation can be expressed. 

Tmin/2 = - (1/2 )Cd^pSv2 

where p = density at given altiude 
V = velocity at minimum drag 
S = wing area of specific airplane 

can be obtained from the above equation, because all values except Cq are 

known, and the drag equation is a function of squared velocity with the known Cr^ . 

9 *^0 

The total value of right side except is the desired constant Kj. 

3. Determining the K 2 

The procedure to compute K 2 is just similiar to that for computing Kj. It is 
convenient to use the same minimum drag concept. At the minimum drag point, the 
minimum thrust required is twice the induced drag amount. 

Use the same maximum endurance velocity and minimum thrust required. 

The induced drag equation is 

^i “ ^min ! ^ 

= (2w2)/(7t e ARp S)(l/v2) 

where w = Gross weight of specific airplane 
s = Wing area of specific airplane 
p = Air density at given altitude 
V = Maximum endurance velocity 

Since the minimum thrust required and all items of the right side are known 
value except l/tt e AR, the value of I/Jt e AR can be solved easily. After obtained the 
l/7t e AR value, substitute this value into the above induced drag equation, and 
calculated the numerical value of the right side except the l/V'^. This numerical value is 
just the desired constant K 2 value. 

But if the raw Kj and K 2 values are used to compute the thrust required, 
there will exist a deviation of more than 10 percent at high speed. Therefore one needs 
to modify the Kj and K 2 to fit the curve to the FLIGHT MANUAL data. 

The fitting method used in this project is the trial and error method. The 
Table 2 and Figure 2.9 illustrates the difference between the modified and raw value of 
Ki,K2- 
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TABLE 2 

THRUST REQUIRED AND FUEL FLOW DEVIATION 
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1 igurc 2.10 Drag Index versus Kl. 



4. Relationship between Kp K 2 and Variation Factors 

Next consideration is to make the equations of K| and K 2 functions of drag 
index, altitude and gross weight. 

Since the K ^ is the function of drag index and air density, K j will be changed 
by different drag index and different altitude. Thus it is needed to develop a 
relationship of K| with the drag index and altitude. 

To accomplish this, it is convenient to select some reference altitude and gross 
weight, because the relationship of K j and drag index is independent of altitude and 
gross weight. Each Kj value for every drag index can be calculated with the previous 
procedure. After calculation of all of the Kj values, a graph of K| versus drag index is 
plotted to find the relationship. 

In this project, the altitude and gross weight are selected to standard sea level 
and 40,000 lbs. The final values and relationships are shown in Table 3 and Figure 
2 . 10 . 



TABLE 3 

DRAG INDEX VS K1 



Drag Index 

0 

20 

40 

60 

80 

100 

120 



0.011534 

0.013335 

0.014807 

0.016780 

0.018480 

0.020527 

0.022361 



Fortunately, as can be seen in the result of the graph, the relationship between 
K| and the drag index is linear. So the equation Kj as a function of drag index is 

Kj = (0.0115 + 0.0001 X Drag Index) 

In this equation, the most important thing is the slope of the graph, because 
the starting value of Kj was modified to fit the published data. 

The Kj equation with modified starting value of Kj is 

Kj = (0.014 + 0,0001 X Drag Index) 

This equation was verified by calculating the drag at each speed. The results 
are shown in Tables 4, 5, 6. 



29 



TABLE 4 

THRUST REQUIRED DEVIATION AT DI 20 



ALTITUDE 


0 ft 


GROSS WEIGHT 40,000 LBS D.I 20 


MILITARY 


THRUST 


21,800 LB 


THRUST 


FUEL FLOW X 0.8489 


FUEL FLOW 


25,680 LB/H 


AIR SPEED 
(KTS) 


FUEL FLOW 
(LBS/H) 


THRUST (LB) 
HAND BOOK 


THRUST (LB) 
CALCULATION 


DEVIATION 

(’^) 


560 


7,725 


6,558 


6,456 


-1.5 


400 


9,014 


7,652 


7,741 


1.2 


440 


10,623 


9,018 


9,198 


2.0 


480 


12,462 


10,579 


10,819 


2.3 


520 


14,652 


12,438 


12,600 


1.3 


560 


17,377 


14,752 


14,536 


-1.5 


600 


21,393 


18,160 


16,625 


-8.5 




FUEL FLOW 


5,600 (LBS/H) 


MAXIMUM I 


ENDURANCE 


THRUST (min) 


4,754 (LB) 




AIR SPEED 


248 (KTAS) 


VARIATION OF K 




0.014+0.0001x20 = 0.016 


*=2 


2.0 X loS 



From Tables 4, 5, 6, the following can be observed. 

First, the thrust required, i.e., calculated drag, is almost same as published 
data at low drag index range, but the thrust required is much less than published data 
at high drag index range. 
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TABLE 5 

THRUST REQUIRED DEVIATION AT DI 60 



ALTITUDE 



FT 



GROSS WEIGHT 



40,000 LBS 



D. I 



60 



MILITARY 


THRUST 


2] 


.,800 lb 


THRUST 


FUEL FLOW X 0.8925 


FUEL FLOW 


24,425 LB/H 


AIR SPEED 
(KTS) 


FUEL FLOW 
(LBS/H) 


THRUST (LB) 
HAND BOOK 


THRUST (LB) 
CALCULATION 


DEVIATION 

C“^) 


360 


9,286 


8,288 


7,935 


- 4.3 


400 


10,980 


9,800 


9,566 


- 2.4 


440 


13,118 


11,708 


11,407 


- 2.6 


480 


15,539 


13,869 


13,449 


- 3.0 


520 


18,662 


16,656 


15,686 


- 5.8 


560. 


23,639 


21,090 


18,115 


1— 1 

• 

r — 1 


600 










MAXIMUM ENDURANCE 

. 


FUEL FLOW 


6,200 (LBS/H) 


THRUST (min) 


5,534 (LB) 


A] 


lr speed 


240 (KTAS) 



VARIATION OF K 



K, 



K, 



0.014+0.0001x60 = 0.02 



2.0 X 10 



8 



Second, in all drag index ranges, the thrust required is much less than 
published data at high air speed range (approximately above Mach number 0.8). 

To compensate for the first problem, that is low thrust required in high drag 
index range, it is needed to increase the Kj value, but the thrust required value is 
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TABLE 6 

THRUST REQUIRED DEVIATION AT DI 80 



ALTITUDE 



FT 



GROSS WEIGHT 



40,000 LBS 



D.I 



80 



MILITARY 


THRUST 


2 


1,800 lb 


THRUST 


FUEL FLOW X 0.9165 


FUEL FLOW 


23,785 lb/H 


AIR SPEED 
CKTS) 


FUEL FLOW 
(LBS/H) 


THRUST (LB) 
HAND BOOK 


THRUST (LB) 
CALCULATION 


DEVIATION 

("O 


360 


10,104 


9,261 


8,674 


- 6.3 


400 


12,094 


11,085 


10,479 


- 5.5 


440 


14,487 


13,278 


12,512 


- 5.8 


480 


17,419 


15,965 


14,763 


- 7.5 


520 


21,477 


19,658 


17,228 


-12.6 


560. 










600 










MAXIMUM ENDURANCE 


FUEL FLOW 


6,650 (LBS/H) 


THRUST (min) 


6,003 (LB) 


A] 


:r speed 


238 (KTAS) 



VARIATION OF K 



K, 



K. 



0.014+0.001x80 = 0.022 



2.0 X 10 



8 



slightly larger than the published data. The way to compensate in both the low and 
high drag index range is to increase the slope of the Kj equation and to decrease the 
starting value. 
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The modified equation of K j is 

Kj = (0.0135 + 0.000132 x Drag Index) 

And the change of Kj equation can be illustrated in Figure 2.11. 

The second problem, that the thrust required is much less than the FLIGHT 
MANUAL data in the speed range higher than Mach number 0.8, can be illustrated as 
the airplane reaches very high flight speeds, the drag rises in a very rapid fashion due 
to compressibility. Since the generalized equation for parasite drag does not account 
for compressibility effects, the actual drag rise is typified by the dashed line in Figure 
2 . 12 . 

This phenomenon can be investigated from the viewpoint of drag coefficient. 
The lift and drag coefficients vary with the Mach number, however, the drag coefficient 
for a fixed value of lift coefficient does not change until the drag-divergence Mach 
number is reached, at which time the drag coefficient increases sharply. If wind-tunnel 
tests of a model of the airplane to be analyzed have been conducted, curves of drag 
coefficient against Mach number for various values of lift coefficient will be available 
as shown in Figure 2.13. 

Drag data above the divergence Mach number cannot be calculated with 
precision, and to achieve any accuracy at all, wind-tunnel tests or other empirical data 
must be relied on. It is noted that even wind-tunnel data may not be reliable in the 
transonic range, and therefore performance guarantees in the transonic range, based 
purely on wind-tunnel data and/or calculations, cannot be considered completely 
dependable. 

In most instances, the lift-drag relation may be expressed in an analytic form 

such as 

Cjjj ~ ^D ^L^ / Jt e AR 

in which case performance characteristics may be analytically determined if similiar 
analytical expression for power available can be written. 

The equation 

Cq ~ ^D ^L^ / e AR 

applies to most aircraft through the Mach number range of 0 to about 4 and a 
range from 0 to 0.6, for both trimmed and untrummed conditions, however, and e 
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Figure 2.11 Drag Index versus Modified KI. 
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Figure 2.12 Typical Airplane Drag Curve. 

become functions of Mach number above Mach numbers of about M = 0.8. For 

supersonic aircraft the value of Cpv increases by a factor of 2 to 3 through the 

0 

transomc region and thereafter generally exhibits a slight decrease. [Ref 3:pp. 275-278] 
The span efficiency factor e shows a steady decrease from M = 0.8 to about M = 2.0 
and thereafter remains almost constant at about 50 to 60% of its low speed value. 
From M = 0.8 to M= 1.2 the parasite drag curve quite closely resembles a sin"^ curve, 
and for approximation analysis may be represented by the following equation 
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DRAG COEFFICIENT IN TRANSONIC 




Figure 2.13 Variation of Cd in Transonic region. 



C 



D 



o 



= Cd + (Cd - Cd ) Sin^ [(M-0.8)/0.8] 

00 om 00 

where Cd = maximum value of parasite drag coefficient 
om 

Cd = minimum low speed value of parasite drag 
00 

coefficient. 
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The variation of e may normally be fitted by the polynomial 

e = eo - e2 ( M - 0.8 + 03 ( M - 0.8 

where e^ = low speed value of e 

^2’ ~ constants which depend on the airplane. 

With the above theory and parasite drag equation, it is possible to make the 
drag equation between the Mach number M.8 to M 1.2. 

Dq = (1/2)pSCd^v2 

at above Mach number M.8, where Mach Drag exist, 

Dq = (l/2)p S [Cq + (Cj) - Cj) )sin3 {(M-0.8)/0.8} ] 

^00 ^om ^00 

Ki = (1/2) pSC^ 

Do = [Ki + Ki / (Cq - Cd ) sin^ {{M-0.8)/0.8} ] 

^ ^00 ^om ^00 

= Kj [ I + (Cp - Cq )/Cd sin^ {(M-0.8)/0.8} ] 

* ^om ^00 ^00 

Since Kj and K2 are function of air density, it must be considered altitude 
correction factor of the changing altitude. [Ref. 4] 

First consideration is the lapse rate and vertical temperature structure for the 
standard atmosphere. 

The adopted value for the lapse rate(-dT/dH = a) in the troposphere is the 
following constant 

a) In terms of standard geopotential meters 

a = 0.0065“C m’^ 

b) In terms of standard geopotential feet (a) 

a = 0.00356616®Ffr^ 

c) In terms of the c.g.s. unit of H 

a = 6.628155 x 10'^ “C cm'^ sec'^ 



Accordingly in the troposphere 



T = Tq - aH 



(b) 



where a,H and Tq are in consistent unit 
Second consideration is the relationship between pressure and vertical 
displacement. 
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The following relationships are to be understood as expressed in any system of 
consistent units. 

The air is assumed to obey the perfect gas law, which, in the c.g.s system of 
units, may be written as 

p = (l/RXP/T) (c) 

where 

p density of air, gm cm'^ 

2 

p pressure, dynes cm' 

T absolute temperature,® K 

R gas constant for 1 gram of dry air, ergs gm'^ {®K}'^ 

And the air is assumed also to be in gydrostatic equilibrium and to satisfy the 
differential equation 

dp = - g dz (d) 

where, in c.g.s units, 

P pressure, dynes cm'^ 
p density (specific mass), gm cm'*’ 

Z vertical distance, cm 
g gravitational acceleration, cm sec*'^ 

The vertical displacement is herein expressed in units of geopotential. 
Geopotential is defined in differential form by the equation 

G dH = g dz (e) 

where 

Z altitude measured positively upwards at a point 
g positive (absolute numerical) value of the acceleration 
due to gravity at the point 
H geopotential at the point 

G dimensional constant, the amount of which determines 
the magnitude of the unit of H in terms of length and 
time.(the dimension of G are in units of gZ per unit 
of H.) 

Substituting equation (b) in equation (c), equation (e) in equation (d), and 
combining the results gives 

dP/P = G/R (-dH/Tg - aH) 
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= G/aR d(To - aH)/(T^ - aH) 



Intergrating the right-hand member between limits of 0 and H gives 

ln( P/Pq ) = G/aR ln( - aH/T^ ) 

Let n = G/aR 
Then 

P/P„ = (T„ - aH/T^ = ( T/Tq (0 

From equation (f), using the consistent numerical values of G,a,and R 
n = 5.2561 (dimensionless) 

The Figure 2.14 shows the relationship between temperature, pressure, density 
and altitude. 

With several equation investigated in above, the useful relationship between 
temperature, pressure change and vertical displacement. 

For this project, the useful relationships are written as 

T/Tq = 1 - 6.875 X 10*^ H 
P/Pq = ( 1 - 6.875 X 10'^ H )5-2561 
T = 518.688 X ( 1 - 6.875 x 10’^ H ) 
a^ = 1.4 X 1714.87 xT 

Finally the relationship between density of air and vertical displacement is 
p/Po = (P/RmPJR^TJ 

= (P/Po)(Tq/T) 

And it can be seen by computing several cases that the effect of changing 
altitude becomes larger than normal density changing ratio. 

Let K 9 =(p/Po)^-2 

Then, Kj is proportional to the K 9 and K 2 is inverse proportional to the K 9 . 

As was shown in previous equations, K 2 is proportional to the squared gross 
weight, and the effect of changing gross weight becomes larger than normal gross 
weight changing ratio, same as altitude change. Thus, if there was a specific reference 
value of K 2 for specific gross weight, the relationship between arbitrary gross weight 
and reference gross weight can be made easily as 
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SPECIFIC VARIATION VS ALTITUDE 




40 



I iljuic 2.14 Icmpciauiic, I’rcssiiic, Density against Gcopotcniial II. 



K2 = K2_.^^(W/W^,f)2-7 

With given three specific conditions, i.e. altitude, gross weight and drag index, 
the thrust required can be computed for any airspeed condition with equations shown 
above. Kj and K 2 values can be computed with specific altitude, gross weight and drag 
index, and the thrust required for any air speed can be computed with K| and K 2 * 

The final goal of this project is to compute the fuel flow rate for some given 
condition, i.e. specific altitude, gross weight, drag index and air speed. The fuel flow is 
thrust required devided by factor C. The factor C is the value that military power 
devided by fuel flow of that time, as was shown in previous step. 

C = MILITARY THRUST / FUEL FLOW (at military thrust) 

5. Factor C 

To compute factor C, it is necessary to find the relationship mihtary thrust 
change for different altitude and fuel flow change. This is done in two steps. 
a. Computing the military thrust 

Altitude is one factor which strongly affects the performance of the turbojet 
engine. An increase in altitude produces a decrease in density and pressure and, if 
below the tropopause, a decrease in temperature. If a typical nonafterbuming turbojet 
engine is operated at a constant RPM and true airspeed, the variation of thrust and 
specific fuel consumption with altitude can be approximated from Figure 2.14. The 
variation of density in the standard atmosphere is shown by the values of density ratio 
at various altitudes. Typical values of the density ratio at specific altitudes are as 
follows 



Altitude, fl; 
Sea level . 

5.000 

10.000 

22,000 

35.000 

40.000 

50.000 



Density ratio 

1.0000 

0.8617 

0.7385 

0.4976 

0.3099 

0.2462 

0.1532 



If the fixed geometry engine is operated at a constant V (TAS) in subsonic 
flight and constant N (RPM) the inlet velocity, inlet ram, and compressor pressure 
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ratio are essentially constant with altitude. An increase in altitude then causes the 
engine air mass flow to decrease in a manner very nearly identical to the altitude 
density ratio. Of course, this decrease in mass flow will produce a significant effect on 
the output thrust of the engine. Actually, the variation of thrust with altitude is not 
quite as severe as the density variation because favorable decreases in temperature 
occur. The decrease in inlet air temperature will provide a relatively greater combustion 
gas energy and allow a greater jet velocity. The increase in jet velocity somewhat 
offsets the decrease in mass flow. Of course, an increase in altitude provides lower 
temperatures below the tropopause. Above the tropopause, no further favorable 
decrease in temperature takes place so a more rapid variation of thurst will take place. 
The approximate variation of thrust with altitude is represented by Figure 2.15 and can 
be computed with the equation as follows 

Tavail = T53,6/.T(l/2) 

where 6 = / P 55 J 

^ ~ Palt ! Pssl 

Some typical values at specific altitudes are as follows 



Altitude, ft; Ratio of (Thrust at altitude/thrust at SSL) 
Sea level 1.0000 

5.000 ,■ 0.8960 

10.000 0.8003 

15.000 0.7114 

20.000 0.6296 

25.000 0.5544 

30.000 0.4855 



Since the change in density with altitude is quite rapid at low altitude, turbojet 
takeoff performance will be greatly affected at high altitude. Also note that the thrust 
at 35,000 ft. is approximately 39 percent of the sea level value. 

The thrust added by the afterburner of a turbojet engine is not affected so 
greatly by altitude as the basic engine thrust. The use of afterburner may provide a 
thrust increase of 50 percent at low altitude or as much as 100 percent at high altitude. 
[Ref 2:pp. 119-121] 
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When the inlet ram and compressor pressure ratio is fixed, the principal 
factor affecting the specific fuel consumption is the inlet air temperature. When the inlet 
air temperature is lowered, a given heat addition can provide relatively greater changes 
in pressure or volume. As a result, a given thrust output requires less fuel flow and the 
specific fuel consumption is reduced. While the effect of altitude on specific fuel 
consumption does not compare with the effect on thrust output, the variation is large 
enough to strongly influence a typical variation of specific fuel consumption with 
altitude. Generally, the specific fuel consumption decreases steadily with altitude until 
the tropopause is reached and the specific fuel consumption at this point is 
approximately 80 percent of the sea level value. 

Above the tropopause the temperature is constant and altitude slightly 
above the tropopause causes no further decrease in specific fuel consumption. Actually, 
altitudes much above the tropopause bring about a general deterioration of overall 
engine efficiency and the specific fuel consumption begins an increase with altitude. 
The extreme altitudes above the tropopause produce low combustion chamber 
pressures, low compressor Reynolds Numbers, low fuel flow, etc, which are not 
conducive to high engine efficiency. 

The available military thrust for this project with various altitude is shown 
on Table 7 and Figure 2.15. 



TABLE 7 

.MILITARY THRUST FOR DIFFERENT ALTITUDES 



Altitude (ft) 

Sea level 

5.000 
l(),000 
15 000 

20.000 
25,000 
30000 



Military thrust (lbs) 

21,800 
19,533 
17,446 
15 508 
13,725 
12,086 
10,584 



The appropriate equation for the curve is as follows 

THRUST = 21,800 - 0.4 x ALTITUDE 

b. Computing the fuel flow at military thrust 

There are three effective factors required to change the fuel flow. One is 
altitude, another is drag index, and the other is gross weight. To get an appropriate 
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MILITARY THRUST VS ALTITUDE 
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I'igurc 2.15 Military 'l lirust against Altitude. 



fuel flow equation as a function of altitude, drag index and gross weight, find the 
relationship between fuel flow and each effective factor, then combine the relationship 
totally. 

(1) Relationship between fuel flow and altitude. 

According to the NATOPS FLIGHT MANUAL data, the fuel flow is 
decreasing gradually with increasing the altitude with other conditions constant. To 
find a relationship between fuel flow and altitude, specific drag index and gross weight 
are fixed, to select the available data. Table 8 and Figure 2.16 show the relationship 
when drag index and gross weight are fixed to zero and 40,000 lbs. 

After curve-fitting, the fuel flow equation stated as follows 

Fj = 26,000 - 0.64 x Altitude (at below 8,000 ft) 

Fj = 20,840 - 0.37 x (Altitude - 8,000 ) (at above 8,000 ft) 
where 

F^ = Fuel flow considered with different altitude 

(2) Relationship between fuel flow and drag index. 

According to the NATOPS FLIGHT MANUAL data, the fuel flow is 
decreased as the drag index increased. As was considered in step (1), a specific altitude 
and gross weight will be fixed to find the relationship between fuel flow and drag index. 

Table 9 and Figure 2.17 show the relationship when altitude and gross 
weight are fixed to zero and 40,000 lbs. 

To find the relationship between fuel flow and drag index with different 
altitude, the computed fuel flow from step (1) must be used. With above two variable 
coefficients, the fuel flow equation is 

^2 “ ^1 * ^ Drag Index (at below Drag Index 20) 

F 2 = (F^ + 333.3) - 22.1 x Drag Index (at above Drag Index 20) 

where 

Fj = Fuel flow came from step (1) 

F 2 = Fuel flow considered with different altitude and 
drag index 

(3) Relationship between fuel flow and gross weight. 

According to the NATOPS FLIGHT MANUAL, the fuel flow 
decreased a small amount with increasing gross weight. To find the relationship 
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FUEL FLOW VS ALTITUDE 
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I'igurc 2.16 Rclalionship bcluecii Fuel F'low and Alliludc. 



TABLE 8 

FUEL FLOW FOR DIFFERENT ALTITUDES 



Altitude (ft) 



Sea level 

4.000 

8.000 
1^000 



Fuel Flow (Ibs/hr) 



26,000 

23,520 

20,840 

19,000 



TABLE 9 

FUEL FLOW FOR DIFFERENT DRAG INDEXS 



Drag Index 

0 

20 

40 

60 

80 

100 

120 



Fuel Flow 

26,000 
25,680 
25 110 
24,425 
23,785 
23,225 
22,745 



TABLE 10 

FUEL FLOW FOR DIFFERENT GROSS WEIGHTS 



Grossweight (lbs) 

35.000 

40.000 

45.000 

50.000 

55.000 

60.000 



Fuel Flow (Ibs/hr) 

26,000 

26,000 

25995 

25,990 

25,985 

25980 



between fuel flow and gross weight, a specific altitude and drag index must be selected 
as in the previous steps. Table 10 and Figure 2.18 show the relationship when altitude 
and drag index are fixed to standard sea level and zero. 

It can be seen from Figure 2.18 that the relationship between fuel flow 
and gross weight is linear. The final equation of fuel flow with different three variables, 
i.e., altitude, drag index and gross weight, can be written as follows 

Fj = F2 -0.002x(Gross weight-40,000) (at heavier than 40,0001bs) 

F3 = F2 (at lighter than 40,000 lbs) 

where 

F2 = Fuel Flow came from step (2) 
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FUEL FLOW VS DRAG INDEX 
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Figure 2.18 Fuel l low versus Gross Weiglil. 



Fj = Fuel Flow considered with dilTerent altitude, 
drag index and gross weight 

Now the military thrust and fuel flow for any altitude, drag index and 
gross weight can be computed with the above step (1) and (2). The desired factor C can 
be computed with the military thrust and fuel flow. 

Factor C = Military thrust/ Fuel Flow 

To compute the fuel flow for any various specific condition is the final 
goal of this project. The fuel flow can be computed with thrust required, i.e. total drag, 
devided by factor C. The total drag and factor C had been computed in previous step 
already. 



FUEL FLOW = THRUST REQUIRED / FACTOR C 

With the above equation, it is possible to find the required fuel flow for 
different specific conditions without interpolating the flight manual data. 
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III. RESULT 



A. BASE LINE CONDITION 

Given condition 

Gross weight 40,000 lbs 

Drag Index 20 

Altitude Standard Sea Level 



1 . Converting velocity versus fuel flow to velocity versus drag 
Military thrust (at sea level) 21,800 lb 
Fuel flow (at military thrust) 25,680 Ib/hr 

Thrust = Fuel Flow x Factor C 
C = Military thrust / Fuel Flow 

at military thrust 

C = 21,800/25,680 = 0.84891 

The thrust required for each velocity can be computed by fuel flow times 
factor C(0. 84891). The Table 11 shows the result. 

TABLE 11 

THRUST REQUIRED FOR EACH SPEED 



VELOCITY(KTAS) 

360 

400 

440 

480 

520 

560 

600 

MIL 



FUEL FLOW(LBS/H) THRUST REQUIRED(LBS) 



7,725 
9,014 
10,623 
12 462 
14652 
17,377 
21 393 
25 680 



6,558 

7,652 

9,018 

10,579 

12,438 

14,752 

18,160 

21,800 



2. Computing the minimum thrust and maximum endurance speed. 

Maximum endurance speed can be obtained from maximum endurance Mach 
number graph in NAPTOPS FLIGHT MANUAL (see Figure 2.7). The maximum 
endurance Mach number is M 0.378. And maximum endurance fuel flow can be 
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obtained from the maximum endurance fuel flow graph in NATOPS MANUAL (see 
Figure 2.8). 

The minimum thrust can be obtained from the maximum endurance fuel flow 
by multiplying with the factor C computed in step 1 . 

maximum endurance fuel flow ; 5,600 Ib/hr 
^x-nin " ^ 0.84891 = 4,753.9 lb 

3. Computing the Kj 

At maximum endurance speed, the thrust required is equal to the total drag, 
that is twice the induced drag or parasite drag, i.e. induced drag and parasite drag are 
the same at this point. 

^min “ ^total “ ^^i “ ^^o 
Dq = Cd (1/2)pv2a 

where Dq = 4,753.9 / 2 = 2,376.9 lb 
p = 0.0023769 lb sec^/ft^ 

V = 0.378 X 1,116.4 ft/sec 
A = 530 ft^ 

= 2Dq / (pV^A) 

- ° = (4,753.9)/(0.0023769 x 422.18^x 530) 

= 0.02117183 

Dq = (0.021 17183)(1/2)(0.0023769)(530)v2 
= 0.01333 

Thus Kj = 0.0133333 

4. Computing the K 2 

Use the same concept as for Kj 

^min ~ ^total ~ ^^i ~ ^^o 
Dq = 2w2/(pAv2)(l/7teAR) 

l/7teAR = 2W2 / (pAV^Dj) 

= (0.0023769)(530)(422.18)^(2376.9)/(2 x 40000^) 

= 0.16677964 
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Dj = (2w2/pA)(l/neAR)(l/v2) 

= (2 X 40,000^ / 0.0023769 x 530)(0.16677964){l/v2) 

Thus 

K2 = 4.2365 X 10^ 

Dtotal “ Thrust required = KjV^ + K 2 /V^ 

These and K 2 values are modified a little to fit the curve of hand book 
data by the trial and error method introduced in Chapter 2. Thus the thrust required 
equation can be written as 

D = 0.01614v2 + 2.0x10^/v2 



But if there were no consideration about the Mach effects, the drag at high 
Mach number, i.e. greater than M 0.8, will be much less than the hand book data. Cj^ 
becomes a function of Mach number above Mach numbers of about M 0.8. In this 
region, the parasite drag curve quite closely resembles a sin curve, and for 
approximation analysis may be represented by the equation 



D - Cq + (Cd - Cq ) sin^ {(M - 0.8)/0.8} 






00 ‘^om 00 ^ 

Thus the drag equation between the Mach number M 0.8 and M 1.2 



Do = 1/2 pSC^y 

Do = 1/2 pS [Cn + (Cd - Cp )sin3{(M - 0.8)/0.8}] 

^ ^00 ^om ^00 

Kj = 1/2 pSCj)^ 

Do = [K, + K,/Cn (Cn - Cd ) sin^fCM - 0.8)/0.8}] 

= Ki[l + (Cq - Cry )/Cn sin2{(M - 0.8)/0.8}] 

^ ^om ^00 ^00 

where Cy^ - Cpv = 0.55 
^om 00 

Cn = 0.018 
^00 

Do = Kj[l + 0.55/0.018 sin^ {(M - 0.8)/0.8}]v2 
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5. Converting to fuel flow from drag 

Fuel flow for each velocity can be computed from thrust required, i.e. total 
drag, devided by factor C. 

C = MILITARY THRUST / FUEL FLOW (at military thrust) 

Military thrust and fuel flow can be computed with the equation that found in 
Chapter 2. 

Military thrust = 21,800 - 0.4 x Altitude 
= 21,800 - 0.4 xO 
= 21,800 

Fuel flow; 

F j = 26,000 - 0.64 x Altitude 
= 26,000 - 0.64 X 0 
= 26,000 

^2 ~ ■ 16.65 X Drag Index 

= 26,000 - 16.65 X 20 
= 25,667 

Fj = F 2 - 0.002 X (Gross weight - 40,000) 

= 26,667 - 0.002 x (40,000 - 40,000) 

= 25,667 

Thus C = 21,800/25,667 
= 0.84934 

Fuel flow = Thrust required / 0.84934 

Table 12 and Figure 3.1 show the results of the calculation and deviation from 
hand book data. 

B. VARIATION NUMBER 1 - GROSS WEIGHT 
Given condition 

Gross weight 50,000 lbs 
Drag Index 20 

Altitude Standard Sea Level 
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THRUST REQUIRED VS VELOCITY 
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320.0 360.0 400.0 440.0 480.0 520.0 560.0 600.0 

VELOCITY (KTAS) BASE LINE CONDITION 



FUEL FLOW VS VELOCITY 
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320.0 360.0 400.0 440.0 480.0 520.0 560.0 600.0 

VELOCITY (KTAS) BASE LINE CONDITION 



The gross weight is changed only from the base line condition. It is necessary to 
consider the change of Kj and K 2 . 



Kj = 1/2 pCp^S 

As you can see from the above equation, K| is independent of gross weight. 
Thus there was no change in Kj. 

But K 2 is a function of gross weight (as discussed in chapter 3). 

K 2 = 2w2/(p S ;t e AR) 



It can be seen that the K 2 is proportional to the squared gross weight, but the 
effects of varing gross weight become more largely. 



Thus K- 



^50,0001b 



= 



^40,0001b 



(50,000/40,000) 



2.7 



8 . n '>^^2.7 



= 2.0x10® X (1.25) 



= 3.653 X 10 



8 



D = 0.0164v 2 + 3.653xl0^/V^ 



Military thrust = 21,800 1 b 
Fuel flow; 

F| = 26,000 - 0.64 X Altitude 
= 26,000 

F 2 = Fj - 16.65 X Drag index 
= 26,000 - 16.65 X 20 
= 25,667 

F 3 = F 2 - 0.002 X (Gross weight - 40,000) 
= 25,667 - 0.002 x (50,000 - 40,000) 

= 25,647 

C = 21,800/25,647 
= 0.85 

Fuel flow = Thrust required / 0.85 



Table 13 and Figure 3.2 show the results of calculation and deviation from hand 
book data. 
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THRUST REQUIRED VS VELOCITY 




fiu (ai) aaainbaa isnani 



60 



Figure 3.2a 1 Itru.sl Required with Gross Weight variation. 



FUEL FLOW VS VELOCITY 
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l iguic 3.2b ITicl Now witii Gross Weight variation. 



C. VARIATION NUMBER 2 - DRAG INDEX 



Given condition 

Gross weight 40,000 lbs 
Drag Index 40 

Altitude Standard Sea Level 



If the drag index was changed from base line condition, what will be effect on the 
result? It is necessary to consider the K j and K 2 equation. 

Kj = 1/2 pSCp^ 

K 2 = 2w 2 / (pStreAR) 



As can be seen, K 2 is a function of gross weight, altitude, wing area and aspect 
ratio, not of drag index, thus K 2 was unchanged. But Kj is a function of altitude, wing 
area and drag index. Thus the appropriate K j value can be computed with equation 



K, = K, x(K, /K, ) 

‘DI20 ‘DI40 ‘DI20 



^DI40 



Now it is necessary to find the relationship of Kj between drag index 20 and 
drag index 40. Furthermore, if the relationship of Kj between any drag index and 
reference drag index, the various Kj for different drag index can be computed. To 
figure out the relationship of Kj for each drag index, let us compute each Kj value. 
The computing method is the same as introduced in the base line condition step. 



D„ - 1/2 pSC[3^v2 
C[3^ = 2D„/(pSv2) 

Kj = 1/2 pSCj) 

For example 
when DI = 0 

Do = Tmin/2 = 4,443.8/2 =2,221.9 

Cj) = 2(2,221.9)/(0.0023769 x530x438.92) 

° =0.018132 

Thus 

K^ = 1/2 pSCp)^ 
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= 1/2 (0.0023769 x 530 x0.018132) 
= 0.011534 



Table 14 and Figure 3.3 show the result and the relationship of Kj for different 
drag index. 



TABLE 14 

K1 VALUE VERSUS DRAG INDEX 



DI 


n 


MAXIMUM ENDURANCE 


THRUST (LB) 




Kl 




FUEL FLOW 


KTAS 


MINIMUM 


Do 


0 


0.8385 


5,300 


260 


4,444 


0.018 


0.012 


20 


0.8489 


5,600 


250 


4,754 


0.021 


0.013 


40 


0.8682 


5,920 


247 


5,140 


0.024 


0.015 


60 


0.8925 


6,250 


241 


5,578 


0.027 


0.017 


80 


0.9165 


6,520 


238 


5,976 


0.029 


0.018 


100 


0.9386 


6,800 


233 


6,383 


0.033 


0.021 


120 


0.9584 


7,050 


230 


6,757 


0.036 


0.022 



The Kj of the base line condition was modified to fit the drag index. As you can 
see in Figure 3.3 and the investigation in Chapter 2, the relationship between Kj and 
drag index is linear and the equation of K j, function of drag index is 

Kj = 0.0135 + 0.000132 x Drag Index 
Thus 

K, = 0.0135 + 0.000132 x 40 
^DI40 

= 0.01878 

Thrust required = 0.01878V^ + 2xl0^/v2 
Military thrust = 21,800 lb 
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K1 VALUE VS DRAG INDEX 
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Figure 3.3 Kl value versus Drag Index. 



Fuel flow; 

Fj = 26,000 

F 2 = (Fj + 333.3) - 22.1 X Drag Index 
= (26,000 + 333.3) - 22.1 x 40 
= 25,449.3 

F 3 = F 2 - 0.002 X (40,000 - 40,000) 

= 25,449.3 

C = 21,800/25,449.3 
= 0.8566 

Fuel flow = Thrust required/0.8566 

Table 15 and Figure 3.4 show the results of calculation and deviation from hand 
book data. 

D. VARIATION NUMBER 3 - ALTITUDE 

Given condition 

Gross weight 40,000 lbs 
Drag Index 20 

Altitude 4,000 ft 

If the altitude will be changed only from the base line condition, how will it affect 
the Kj and K 2 ? 

The equations of K j and K 2 are 
Kj = 1/2 pSCj3^ 

K2 = 2w2/(pS7teAR) 

As can be seen in the above equation, the Kj and K 2 both are functions of the 
air density, and air density is a function of altitude. Thus the Kj and K 2 must be 
changed with varying altitude. 

The useful relationships came from the ICAO Report written as 
p = P/RT 

T/Tq = 1 - 6.875xl0'^H 
P/Pq = (1 - 6.875x10'%)^-2^^^ 
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TABLE 15 

THRUST REQUIRED AND FUEL FLOW WITH DILAG INDEX VARIATION 
® FT I GROSS WEIGHT I 40,000 LBS I DRAG INDEX 
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THRUST REQUIRED VS VELOCITY 
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320.0 360.0 400.0 440.0 480.0 620.0 560.0 

VELOCITY (KTAS) DRAG INDEX VARIATION 



FUEL FLOW VS VELOCITY 
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320.0 360.0 400.0 440.0 480.0 520.0 560.0 

VELOCITY (KTAS) DRAG INDEX VARIATION 



T = 518.688 x(l - 6.875xlO'^H) 

And the relationship between density of air and vertical displacement is 

P/Po = (P/RTXPq/RTo) 

= (P/PoKT^/T) 

Since the effects of changing altitude become larger, it is necessary to modify the 
equation. 

Let Kg = (p/Po)l-2 

Then is proportional to the Kg and K 2 is inverse proportional to the Kg. 

Thus in this case 

P/Po = (P/PoXTo/T) 

= ( 1 - 6.875xl0’^x4, 000)^-2561/(1 . 6,875x10*5 x 4,000) 

= 0.88809 

K, = K, X (0.88809) P2 

1 1 

= .0139976 

K 2 = K 2 j / (0.88809)P2 

= 2.306x10^ 

Thrust required = 0.0139976v2 + 2.306xl0^/v2 
The military thrust will be reduced as discussed in Chapter 2. 

Tavail = 8/a<'/2) 

^ = P 4 , 000ft ! Pssl 

= 1,827.69/2,116.22 
= 0.8636578 

^ "" P 4 , 000ft / Pssl 

= 0.0021109/0.0023769 
= 0.8880895 

Tavail = 21,800 x0.8636578/0.8880895^/2 

= 19,779 



69 



Fuel flow; 

Fj = 26,000 - 0.64 x Altitude 
= 26,000 - 0.64 X 4,000 
= 23,440 

F 2 = Fj - 16.65 X Drag Index 
= 23,440 - 16.65 x 20 
= 23,107 

Fj = F 2 - 0.002 X (Gross weight - 40,000) 

= 23,107 - 0.002 X (40,000 - 40,000) 

= 0.86463 

C = 19,979/23,107 
= 0.86463 

Fuel flow = Thrust required/0.86463 

Table 16 and Figure 3.5 show the results of calculation and deviation from hand 
book data. 



E. VARIATION NUMBER 4 - WEIGHT, DRAG INDEX AND ALTITUDE 
Given condition 

Gross weight 50,000 lbs 
Drag Index 40 

Altitude 8,000 ft 



If a'l conditions are changed from base line conditions, it is neceassary to 
consider for each changed condition step by step, same as previous sections. 

At first, think about gross weight change. The changing gross weight affects the 
k 2 only, because K j is independent of gross weight. 

K 2 = 2w 2 / (pSJteAR) 



As you can see in section 3.2 






^50,0001bs 



= 



^40,0001bs 
= 2.0x10^ (1.25)2-'7 
= 3.653x10^ 



(50,000/40,000)2-'^ 
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THRUST REQUIRED VS VELOCITY 
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320.0 360.0 400.0 440.0 480.0 520.0 560.0 600.0 

VELOCITY (KTAS) ALTITUDE VARIATION 



FUEL FLOW VS VELOCITY 
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I'iguic 3.5b ITicI l low with Altitude variation. 



Next, think about the drag index change. If the drag index is changed, it is 
necessary to consider K j only, because the K 2 is independent of the drag index change. 
K j equation, function of drag index is 

K| = 0.0135 + 0.000132 x Drag Index 

Thus 

Ki = 0.0135 + 0.000132 x 40 
‘DI40 

= 0.01878 

And finally, think about the altitude change. 

As you can see in section 3.4, the altitude change affects both Kj and K 2 , 
because Kj and K 2 are function of air density. The equation in section 3.4 can be used 
in this step. 



P/Pq = (P/RTKPq/RTq) 

= (P/Po)(Tq/T) 

Let K 9 = (p/Po)‘-2 
Kj = f(p) 

K 2 = f (1/p) 

Thus 

^1 " ^lDi40^^P8,000ft/Pssl)^'^ 

= 0.0 1 878 X (0.00 1 8683/0.0023769) ^ 

= 0.01406755 

~ ^^50,0001bs ^ ^^ 8 , 000 ft/Pssl) 

= 3.653x10^ / (0.0018683/0.0023769)^-^ 

= 2.9157x10^ 

Thrust required = 0.01407V^ + 2.9157xl0^/V^ 

The military thrust will be computed as discussed in Chapter 2. 

^ 8 , 000 ft " ^ssl 

^ = ^ 8 , 000 ft / ^ssl 
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= 1571.88/2116.22 
= 0.742777 

^ ^ P8, 000ft / Pssl 

-0.0018683/0.0023769 
= 0.786024 

”^8, 000ft "" 21,800 X (0.742777/0.786024 1 '2) 

= 18,2641b 

Fuel flow; 

Fj = 26,000 - 0.64 X Altitude 
= 26,000 - 0.64x8,000 
= 20,880 

F 2 = (Fj + 333.3) - 22. 1 X Drag Index 
= 21,213.3 - 22.1 x40 
= 20,329.3 

Fj = F 2 - 0.002 X (Gross weight - 40,000) 

= 20,329.3 - 0.002 x (50,000 - 40,000) 

= 20,309.3 

C = Military thrust / Fuel flow 
= 18,264/20,309.3 
= 0.8993 

Fuel flow = Thrust required/0.8993 

Table 17 and Figure 3.6 show the results of calculation and deviation from hand 
book data. 

F. VARIATION NUMBER 5 - WEIGHT, DRAG INDEX AND ALTITUDE 

We have seen how the fuel flow was changed with various conditions. Not much 
time is needed to fmd the fuel flow for conditions which can be read directly from hand 
book, but most of the conditions do not correspond to the published conditions which 
can be read directly from the hand book data sheets. Thus to find the fuel flow for 
specific conditions that are not printed in the hand book, one must spend lots of time 
to interpolate each item. The specific fuel flow will now be computed by both methods. 
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TABLE 17 

THRUST REQUIRED AND FUEI. I'LOW WITH ALL CONDITIONS 

VARIATION 
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THRUST R1-:QU1RR!) vs VRLOCITY 
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Hguic 3.6a lluiisi. Rc(|uiic(I with All C\)iKlition variation. 



FUEL FLOW VS VELOCITY 
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440.0 480,0 520.0 560.0 

VELOCITY (KTAS) ALL CONDITION VARIATION 



Given condition 

Gross weight 48,000 lbs 
Drag Index 27 

Altitude 1 1 ,000 ft 

Velocity 455 KTS 

The procedure to compute the Kj and K 2 is the same as in the previous section. 
Consider each different condition. 



1. Gross weight 



K- 



M8,0001bs 



= K- 



-40,0001bs 



X (48,000/40,000) 



2.7 



= 2 . 0 x 10 ^ x(1.2)2-7 



= 3.27x10 



8 



K j is independent of changing gross weight. 
2. Drag index 

Kj = 0.0135 + 0.000132 x Drag Index 
= 0.0135 + 0.000132 x 27 
= 0.01706 



K 2 is independent of changing Drag Index. 

3. Altitude 

= ^ 1 qj 27 ^^*^ll»000ft / Pssl)^'^ 

= 0.01706 X (0.0017007/0.0023769)^-2 

= 0.01142 

*^ 2 " K2_,joooib'‘'’ll.O(X)ft/Pssl)‘-^ 

= 3.27x10^/(0.0017708/0.0023769)^-2 

= 4.8893x10^ 

Thrust required = 0.01142v2 + 4.8893xl0^/v2 

= 0.01142 X (455x1.69)2 + 4.8893x10^/(455x1.69)2 
= 7,572 lb 



Compute the factor C; 
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Military thrust = 21,800 - 0.4 x Altitude 
= 21,800 - 0.4x11,000 
= 17,400 1b 

Fuel flow; 

F| = 20,840 - 0.37 x (Altitude - 8,000) 

= 20,840 - 0.37 X (1 1,000 - 8,000) 

= 19,730 

F 2 = (Fi +333.3) - 22.1 X Drag Index 
= 20,063.3 - 22.1 X 27 
= 19,466.6 

F 3 = F 2 - 0.002 X (Gross weight - 40,000) 
= 19,466.6 - 0.002 x (48,000 - 40,000) 

= 19,450.6 

C = Military thrust / Fuel flow 
= 17,400/19450.6 
= 0.89457 



Thus 

Fuel flow = Thrust required / C 
= 7,572/0.89457 
= 8,464 Ibs/hr 

Now the fuel flow for given conditions will be computed by interpolating the 
hand book data. Figure 3.7 shows the whole diagram to compute the desired fuel flow 
of given conditions. 

Computations for the steps are: 

1) FI = (8873 -7988)x(27 - 20)/(40 - 20) + 7988 = 8297.75 

2) F2 = (10324 - 9223)x(27 - 20)/(40 - 20) + 9223 = 9608.35 

3) F3 = (455 - 440)x(9608.35 - 8297.75)/(480 - 440) + 8297.75 = 8789.23 

4) G1 = (7729 - 6962)x(27 - 20)/(40 - 20) + 6962 = 7230.45 

5) G2 = (8840 - 7933)x(27 -20)/(40 -20) + 7933 = 8250.45 

6 ) G3 = (455 - 440)x(8250.45 - 7230.45)/(480 - 440) + 7230.75 = 7613.25 

7) L4 = (11000 - 8000)x(7613.25 - 8789.23)/( 12000 - 8000) + 8789.23 = 7907.24 

8 ) HI = (9239 - 8366)x(27 - 20)/(40 - 20) + 8366 = 8671.55 
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DIAGRAM FOR DESIRED FUEL FLOW 



GROSS WEIGHT 40,000 LBS 


ALTITUDE 8,000 FT 


11,000 FT 


ALTITUDE 12,000 FT 


A/S 

KTAS 


DRAG INDEX 


^4 


A/S 

KTAS 


DRAG INDEX 


20 


27 


40 


20 


27 


40 


440 


7988 


Fl 


8873 


440 


6962 


Gl 


7729 


455 




^3 




455 




^3 




480 


9223 


CN 


10324 


480 


7933 




G2 


8840 









GROSS 


WEIGHT 50,000 LBS 






ALTITUDE 8,000 FT 


11,000 FT 


ALTITUDE 12,000 


FT 


A/S 


DRAG INDEX 




A/S 


DRAG INDEX 


KTAS 


20 


27 


40 




KTAS 


20 


27 


40 


440 


8366 


Hi 


9239 




440 


7447 


c — i 


8192 


455 




H3 




^4 


455 




J3 




480 


9522 


«2 


10688 




480 


8315 


J2 


9216 



G/WT 


40,000 LBS 


48,000 LBS 


G/WT 


50,000 LBS 


A/S 


DRAG INDEX 




A/S 


DRAG INDEX 


KTAS 


27 




KTAS 


27 


455 


^4 


^4 


455 





Figure 3.7 Diagram for Desired Fuel Flow. 
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9) H? = (10688 - 9522)x(27 - 20)/(40 -20) + 9522 = 9930.1 

10) H3 = (455 - 440)x(9930.1 - 8671.55)/(480 -440) -I- 8671.55 = 9143.51 

11) J1 = (8192 - 7447)x(27 - 20)/(40 - 20) + 7447 = 7707.75 

12) J2 = (9216 - 8315)x(27 - 20)/(40 - 20) + 8315 = 8630.35 

13) J3 = (455 - 440)x(8630.35 - 7707.75)/(480 - 440) -I- 7707.75 = 8053.73 

14) • 8000)x(8053.73 - 9143.51)/(12000 - 8000) + 9143.51 

15) N4^ 39 ■ 40000)x(8326.17 - 7907.24)/(50000- 40000) -I- 7907.24 



Desired fuel flow for given condition is 8,242 Ibs/hr. 
Thus deviation of the results 

(8,464 - 8,242)/8,242 x 100 = 2.7 % 
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IV. CONCLUSION 



It is concluded from the investigations that computing the fuel flow rate from the 
computer program results in less deviations and is simpler than interpolating from 
NATOPS FLIGHT MANUAL DATA. The total lines used in this program is 53 lines 
for the m;.in frame program and 229 lines for the hand-held calculator program, and 
both programs are inserted in Appendix. 

Though the deviation of fuel flow for this program was restricted to a maximum 
of 10 percent, for most conditions deviations are distributed below 5 percent. To 
determine the fuel flow at specific conditions which can not be read directly from the 
flight manual, it will be necessary to calculate fifteen steps of interpolations. This is 
very time consuming and the probability of making a mistake is very high when 
calculating so many interpolations. 

Thus the computer program to compute the fuel flow rate can not only save time 
for many conditions which must be changed continuously during flight, but it also 
reduces the probability of making a mistake in interpolations steps. 

It is recommended that a similar computer program be developed for specific 
flight stages, i.e., take-off, climb, or landing, using the same concepts that were used in 
this program. 
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APPENDIX A 

COMPUTER PROGRAM (BASIC) 



00010 K2=2.0E+08 ! REM XX BASED ON S.S.L. AND GROSS WEIGHT AOOOOLBS XX 
00020 PRINT 'INPUT ALTITUDE, GROSS WEIGHT, DRAG INDEX, AIR SPEED, A/S INC' 
00030 INPUT A, W, 19, VI, VP 

00040 Kl=(0.0135+. 000132X19) ! REM XX K1 VARIES WITH DRAG INDEX xx 
00050 D4=l-6 .875E-6XA ! REM xx ICAO SCALE BASE xx 

00060 T1=518.688XD4 ! REM XX ABSOLUTE TEMP (R) BASED ON ALTITUDE XX 
00070 D5=D4XX5.2561 ! REM xx PRESSURE RATIO BASED ON ALTITUDE xx 
00080 A9=(I .4X1714 .87XT1) ! REM xx SPEED OF SOUND (SQUARED) xx 
00090 K9=(D5/D4)XX1.2 ! REM XX ALTITUDE CORRECTION FOR 'Q' XX 
00100 K1=K1XK9 ! REM XX PAXMXX2 ALTITUDE CORRECTION XX 
00110 K2=K2/K9 ! REM XX PAxMxx2 ALTITUDE CORRECTION xx 

00130 K2=K2X(W/40000)xx2.7 !REMXx WEIGHT OTHER THAN 40000LB CORRECTIONXX 
nm/^n print = *k7=**K7 

00170 PRINT TABd); 'VEL(KTS)' ;TAB(12); 'MACH' ;TAB(20); 'PAR DRAG'; 

00180 PRINT TAB(34);'PRO DRAG' ;TAB(48 ); 'MACH DRAG' ; TAB(62 );' FUEL FLOW' 
00190 V2=VIX1 .68894 
00200 M=V2/A9XX0.5 

00210 D8=l : REM XX CDO MACH CORRECTION BASE XX 
00220 IF M>.8 THEN 480 

00230 D1=K1XV2XX2 ! REM XX PARASITE DRAG (SUBSONIC) XX 
00240 D3=D1XD8 ! REM xx PARASITE DRAG (MACH CDO) xx 
00250 D2=K2/V2XX2 ! REM XX INDUCED DRAG XX 
00260 D0=D3+D2 ! REM xx MACH DRAG xx 
00270 D7=DI+D2 ! REM XX SUBSONIC TOTAL DRAG xx 

00280 T=21800-AX0.4 ! REM XX MIL THRUST OTHER THAN S.S.L. CORRECTION XX 
00290 IF A<=8000 THEN 320 

00300 FI = 20840-(A-8000)X0.37 ! REM xx F/FLOW CORRECTTION ABOVE 8000 FTXX 
00310 GO TO 330 

00320 F1=26000-0.64000XA ! REM XX F/FLOW CORRECTION BELOW 8000 FT XX 
00330 IF I9>20 THEN 360 

00340 F2=F1-16 .65X19 ! REM XX F/FLOW CORRECTION BELOW DRAG INDEX 20 XX 
00350 GO TO 370 

00360 F2=(Fl+333.3)-22. 10X19 ! X F/FLOW CORRECTION ABOVE DRAG INDEX 20 X 
00370 IF W<=40000THEN 400 

00380 F3=F2-(W-40000)X0 .002 ! REM XX F/FLOW CORRECTION WITH WEIGHT XX 
00390 GO TO 410 
00400 F3=F2 

00410 C=T/F3 ! REM xx COEFFICIENT FOR FUEL FLOW XX 
00420 FF=D0/C ! REM xx FUEL FLOW XX 

00430 PRINT TAB(3);V1;TAB(8);M;TAB(18);D2;TAB(32);D3;TAB(46);D0; 

00440 PRINT TAB(60);FF 
00450 V1=V1+VP 
00460 IF Vl>610 THEN 510 
00470 GOTO 190 

00480 D9=.55x(SIN((M-.8)/.8)xx3)/.0180 

00490 D8=l+D9 

00500 GOTO 230 

00510 PRINT 'C=';C, 'T=';T 

00520 STOP 

00530 END 
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APPENDIX B 

COMPUTER PROGRAM (HAND-HELD CALCULATOR) 



f'i*LEL ‘FF‘ 


51 RCL 03 


101 RCL 15 


15! + 


201 PCL 23 


e: 1 


52 » 


102 RCL 13 


152 RCL 19 


202 -FF=- 


E* .^70 £'4 


53 STO 09 


103 + 


153 + 


203 ORCL X 


F4 PCL ee 


54 RlL 04 


104 STO 17 


154 STO 20 


204 OVIEH 


f-5 


55 RCL 09 


105 RCL 00 


155*LBL -EE- 


205 STOP 


H6 PRORPT 


56 ' 


106 .4 


156 RCL 01 


206 -HEH V’<1,0> 


y? STO 00 


57 STO 04 


107 * 


157 40000 


207 PROMPT 


00 RCL 01 


5b 2.0 E03 


leS CHS 


15S - 


203 X=8? 




59 RCL 09 


109 21800 


159 .002 


209 STOP 


10 PRORPT 


60 / 


110 + 


160 * 


210 GTO -VV 


11 STO 01 


61 STO 10 


111 STO 13 


161 CHS 


211 ♦LBL -HN- 


12 RCL 02 


62 RCL 01 


112 8000 


162 RCL 20 


212 RCL 12 


13 •IH=?’ 


63 40000 


113 RCL 00 


163 + 


213 .8 


14 PRORPT 


64 / 


114 X<=Y? 


164 STO 21 


214 - 


15 STO 02 


65 2,7 


115 GTO -PP- 


165 1/X 


215 .8 


It^LBL -VV 


66 YtX 


116 RCL 00’- 


166 PCL 18 


216 / 


17 RCL 03 


67 RCL 10 


117 8000 


167 * 


217 EHTERt 


IS •V<KT3>=?- 


6S » 


118 - 


163 STO 22 


218 3 


19 PPORPT 


69 STO 10 


119 .37 


169 1/X 


219 YTX 


20 STO 03 


70 RCL 03 


120 » 


170 RCL 16 


220 SIN 


21 RCL 02 


71 1.68394 


121 CHS 


171 » 


221 .55 


22 .000132 


72 * 


122 20840 


172 STO 23 


222 * 


23 » 


73 STO 11 


123 + 


173 PCL 03 


223 .018 


24 .0135 


74 RCL 07 


124 STO 19 


174 FIX 0 


224 / 


25 + 


75 SORT 


125 GTO -LL- 


175 •va;TS>= 


225 1 


20 STO 04 


76 I/X 


126 ♦LBL -PP- 


176 CRCL X 


226 + 


27-F..S75 E-00 


77 RCL 11 


127 PCL 00 


177 OVIEH 


227 STO 24 


23 RCL 00 


78 * 


12E .64 


178 STOP 


223 GTO -FF- 


29 » 


79 STO 12 


129 * 


179 RCL 12 


229 .EHD. 


30 CHS 


80 .8 


130 CHS 


160 FIX 3 




31 1 


81 X^=Y? 


131 26000 


181 -MCH=- 




32 + 


82 GTO -HR* 


132 + 


182 PRCL X 




33 STO 05 


33»LeL -RR- 


133 STO 19 


183 PVIEH 




34 518.08 


04 RCL 11 


134»L6L -LL* 


184 STOP 




35 * 


85 Xt2 


135 26 


185 RCL 15 




36 STO 06 


86 RCL 04 


136 RCL 02 


186 FIX 1 




37 1716.5 


87 * 


137 X)Y? 


187 -D<I)=- 




33 » 


88 STO 13 


138 GTO -PD- 


183 flRCL X 




39 1.4 


89 RCL 24 


139 16.65 


189 OVIEH 




40 ♦ 


90 * 


140 ♦ 


190 STOP 




41 STO 07 


91 STO 14 


141 CHS 


191 RCL H 




42 RCL 05 


92 RCL 11 


142 RCL 19 


192 *D<P>=- 




43 EHTEPt 


93 Xt2 


143 ♦ 


193 ORCL X 




44 5.2561 


94 1/X 


144 STO 20 


194 OVIEH 




45 YtX 


95 RCL 10 


145 GTO -EE- 


195 STOP 




46 STO 03 


96 * 


146*LBL -DD’ 


196 RCL 16 




47 1116 


97 STO 15 


147 22.1 


197 -D<f1>=* 




43 Xt2 


98 RCL 14 


148 ♦ 


198 ORCL X 




49 RCL 07 


99 ♦ 


149 CHS 


199 OVIEH 




50 / 


100 STO 16 


150 333.3 


200 STOP 





85 



LIST OF REFERENCES 



1. Kuethe, Arnold M. and Chow. Chuen-Yen Foundation Of Aerodynamics, 3rd 

ed., JoHn Wiley & Sons, Inc., 1^76. 



2. Hurt, H. H. Jr., Aerodynamic For Naval Aviators, OfTice of the Chief of 

Naval Operations Aviation Training Division, 1965. 



3. Domash, Doniel O. Sherby, Sydney S. and Connolly, Thomas F. Airplane 

Aerodynamics, 3rd ed.. Pitman, 1961. 



4. NASA Report 1235, National Advisory Committee For Aeronautics, 1955. 



86 



INITIAL DISTRIBUTION LIST 



No. Copie.s 



Defense Technical Information Center 

Carmeron Station 

Alexandria, Virginia 22304-6145 

Library, Code 0142 
Naval 'Postgraduate School 
Monterey, California 93943-5002 

Department Chairman, Code 67 
Dept, of Aeronautics 
Naval Postgraduate School 
Monterey, California 93943-5000 

Prof Donald M. Layton, Code 67Ln 
Dept, of Aeronautins 
Naval Postetaduate School 
Monterey, California 93943-5000 

Oh, Chang Whan 

Yun-Hee 2 Dong 95-14, Seo dae Moon-ku, 
Seoul, Republic of Korea 

Air Force Central Library 
Sindaebang Dong. Kvvanak Gu, 

Seoul, Republic of Korea 

3rd Department of Air Force College 
Sindaebang Dong. Kwanak Gu, 

Seoul, Republic of Korea 

Library of Air Force Academy 
Chongwon Gun, Chung Cheong Bug Do, 
Republic of Korea 



2 



2 



1 



1 



10 



2 



2 



2 



r 



87 



18070 2 




^ ' '' W "I ‘ij^lM :yIB ..; 










2211i,0 ’ 

Thesis 
034275 Oh 

c.l The use of a computer 

to obtain flight manual 
data. 

2 OCT C9 ^0 0, Cl 6 0 > 

* J f LB O ^ 



221140 

Thesis 
034275 Oh 

c.l The use of a computer 

to obtain flight manual 
data. 



